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1.  Introduction 


Let  G  be  a  finite  dimensional  Lie  group.  In  this  paper,  we  consider  left-invariant 
control  systems  on  G.  To  fix  notation,  let  Lg  :  G  — »•  G,  h  t— »  g  ■  h,  denote  the  left 
translation  by  g  £  G,  and  let  ThLg  :  ThG  t— >  TghG  denote  its  tangent  (i.e  linearization). 
Let  Q  denote  the  Lie  algebra  of  G  (i.e  TeG  the  tangent  space  at  identity  e  €  G).  A  vector 
field  g  vg  €  TG  is  left-invariant  if  vg  =  TeLg  •  £  for  some  fixed  £  £  Q ■  A  left- invariant 
control  system  is  defined  by  letting  £  be  a  (controlled)  curve  in  Q.  It  takes  the  form 


g  =  TeLg  •  £u  (1) 

where  each  control  «(•)  determines  a  curve  £u(.)  C  Q.  Here  we  limit  ourselves  to  vector 
valued  control  functions  «(•). 


=  £o  +  XI Ui  ^  (2) 

i=i 

where  Uj  ( t )  €  H,  {£o,  £i,  •••,  £m}  spans  an  ( m  +  1  )-dimensional  subspace  7i  of  £?, 
m  +  1  <  n  —  dim  (Q). 

Consider  an  optimal  control  problem  of  the  form 

T 

Min  f  .  , 

/  L(u)  d.t  (3) 

«(•)  J 

o 

subject  to  the  condition  that  u(-)  steers  (1)  -  (2)  from  go  at  t  —  0  to  g\  at  t  =  T.  (T 
may  be  fixed  or  free.  For  simplicity  we  restrict  T  to  be  fixed  in  this  paper.) 

Clearly,  the  Lagrangian  L  is  G-invariant  (-does  not  depend  on  g). 

It  is  the  content  of  the  maximum  principle  that  optimal  curves  in  G  are  base  integral 
curves  of  a  hamiltonian  vector  field  on  T*G.  To  be  more  precise,  let  tq  :  TG  — *■  G  and 
Tq  :  T*G  ->  G  be  bundle  projections. 
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Define 


Hx  =  Hx  ( ag ,  u) 

=  -  A L(u)  +  <  ag,  TeLg  ■  > 


(4) 


where  A  =  1  or  0,  and  ag  £  T*G. 

Maximum  Principle 

Let  uopt  be  a  minimizer  of  (3)  and  let  g  be  the  corresponding  trajectory  of  (1).  Then, 
g(t)  =  Tq  (ag(t))  for  an  integral  curve  ag  of  the  hamiltonian  vector  field  X^lpt  defined 
for  t  £  [0,  T]  such  that: 

(a)  If  A  =  0  then  ag  is  not  the  zero  section  of  T*G  on  [0,  T]. 

(b)  HA(a9,  uopt )  =  suPjHx(ag,  u )  for  t  a.e.  in  [0,  T}.  Here  U  is  the  space  of  values  of 
controls  =  IR,”'  in  this  paper. 

(c)  If  the  terminal  time  T  is  fixed  then  HA(a9,  uopt)  —  constant,  and  if  T  is  free  then 
Ha(gi9,  uopt )  =  0  V<  £  [0,  T).  Trajectories  corresponding  to  A  =  0  are  called  abnormal 
extremals  and  they  occur  often  but  are  ruled  out  under  suitable  hypotheses.  We  are 
concerned  solely  with  regular  extremals  (\  =  1)  in  this  paper.  In  that  case,  since  U  —  ]Rm, 
optimal  controls  for  (3)  subject  to  (1)  and  (2)  satisfy 

1h~  +  d^i  {<X 9'  TeL9'iu)  =°-  (5) 

*  =  1,2, 


From  (2), 


(otgi  TeLg  •  £ u )  —  (agi  TeLg  ■  (£0  +  'y  ^  uiG)) 

i=  1 
m 

=  ( TeL*-ag ,  £0  + 

2=1 

m 

=  (/b  Co)  +  it)  (6) 

i= 1 
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where  p,  —  TcL*  •  ag  €  Q* .  From  (5)  and  (6) 


-dL 

dui 


+  (/b  £i) 


i 


0, 

1, 2,  •  •  • ,  m. 


7) 


r  A=1 


are 


Observe  that,  from  (6)  and  (7)  the  pre-hamiltonian  Hx  1  and  the  hamiltonian  H 

rn 

G-invariant.  More  explicitly,  suppose  L(u)  =  ^  i*  u],  the  constants  U  >  0.  Then,  the 


optimal  Ui  are  given  by 


i=l 


and  the  hamiltonian  on  T*G  is 


Ui 


(ft>  Zi) 

Ii 


(8) 


H  =  5.)  +  \  ±  ^  (9) 

1=1 

clearly  the  hamiltonian  H  is  G-invariant.  This  brings  us  to  the  main  (and  essentially 
only  •point )  of  this  paper.  The  hamiltonian  vector  field  Xh  on  T*G  corresponding  to 
the  hamiltonian  H  can  be  (left-)  reduced  to  a  hamiltonian  vector  field  Xh  on  Q* .  The 
latter  is  hamiltonian  in  a  non-canonical  (Lie-Poisson)  sense.  Integral  curves  of  X  h  can 
be  reconstructed  from  integral  curves  of  Xj,  by  quadrature.  Thus  questions  about  explicit 
solvability  (or  integrability)  of  Xh  are  turned  into  corresponding  questions  about  Xh- 

The  process  by  which  one  passes  from  Xh  to  Xh,  known  as  Poisson  reduction,  is  well- 
known  in  geometric  mechanics.  Since  Q*  is  a  vector  space,  we  thus  conveniently  bypass 
computations  with  the  symplectic  structure  on  T*G.  In  his  papers  Jurdjevic  works  with 
the  symplectic  structure  on  T*G  although  one  does  see  reduction  “done  by  hand”  in  special 
examples. 

The  paper  of  Faybusovich  is  closest  to  the  point-of-view  of  the  present  paper  and  is  in 
some  ways  broader  in  scope.  Limiting  ourselves  to  G-invariant  optimal  control  problems 
as  here,  we  hope  to  make  things  a  bit  more  explicit  than  is  perhaps  in  [2]. 
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In  the  next  few  sections  we  collect  together  some  basic  facts  regarding  Poisson  man¬ 
ifolds,  summetry  and  reduction  and  carry  out  the  reduction  step  for  the  optimal  control 
problem  (l)-(3). 

2.  Poisson  Manifolds,  Symmetry  Groups  and  Reduction 

A  smooth  manifold  P  is  said  to  be  a  Poisson  manifold  if  it  comes  equipped  with  an 
IR-bilinear  map  (Poisson  structure), 


{-,-}p  :  C°°(P)  x  C°°(P)  ->  C°°(P) 

satisfying  the  axioms: 

(*){^»  V> }p  =  -  {Vb  <i>) 

(«W>  {Vb  x}p}p  +  {Vb  {x,  </> }p}p  +  {x,  {<f>,  4'}p}p  =0 
(iii){H>,  x}p  =  *l>  {<!>,  x}p  +  tHVb  X } p 

Associated  to  a  Poisson  structure  is  a  unique  twice-  contravariant,  skew-  symmetric, 
tensor  field  A  on  P  on  such  that 


{<j>,  %1)}P  =  A (d<f>,  d*l>).  (10) 

For  a  proof  see  p.  109  of  [8].  The  tensor  A  defines  a  vector-bundle  morphism, 

A#  ;  t*P  -►  TP 

ax  A*(ax)  €  TXP 

satisfying, 

f3x  (A#  (ax))  =  A(x)  (ax,  /3X)  for  all  ax,  f3x  £  TXP. 

The  rank  of  the  Poisson  structure  at  x  £  P  is  defined  to  be  the  rank  of  the  Poisson  ten¬ 
sor  A  at  x.  This  is  simply  the  rank  of  the  (characteristic)  distribution  C  —  A #(T*P)  C  TP 
at  the  point  x.  The  rank  may  vary  on  P.  It  is  a  theorem  of  Kirillov  [6]  that  A#  (T* P) 
defines  a  generalized  foliation  on  P  such  that  through  each  point  x  £  P,  passes  a  leaf 
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carrying  a  unique  symplectic  structure  that  makes  the  injection  map  of  the  leaf  a  Poisson 
morphism.  (See  Weinstein  [13]  and  Liebermann-Marle  [8]).  Thus  a  Poisson  manifold  is  a 
union  of  symplectic  leaves. 

A  function  <j>  G  C°°(P)  is  called  a  Casimir  function  if  {0,  ift}p  =  0  V?/’  G  C°°(P). 
Casimir  functions  are  constant  on  symplectic  leaves. 

Let  G  be  a  Lie  group  and  let  :  G  x  P  — >  P,  (g,  x)  i— >  Ts(x)  be  a  group  action 
such  that,  ’!'(,(•)  is  a  Poisson  morphism  for  every  g  G  G.  Further  suppose  that  the  action 
is  proper  and  free.  Then  the  quotient  manifold  P / G  is  smooth  and  carries  the  Poisson 
structure  {•,  -}p/g  induced  from  the  one  on  P  satisfying. 


{<j>,  fp}p/G°^  =  0  ott}p, 

where  7 r  :  P  P/G  is  the  canonical  projection.  By  construction,  it  is  a  Poisson  morphism. 

Now  G-invariant  dynamics  on  P  induces  dynamics  on  P/G.  To  see  this,  let  H  :  P  — » 
It  be  a  G-invariant,  hamiltonian  function  on  P,  i.e, 

H(*g(x))  =  H(x)  VyeG. 

Define  the  (hamiltonian)  vector  field  Xp  on  P  by  requiring, 


XH  f  =  {<f>,  H}  M  f  G  C°°(P). 

The  hamiltonian  H  descends  to  h  :  P/G  — »  1R  and  determines  the  reduced  hamiltonian 
dynamics  X T  on  P/G  by 


Xhf=  {l  h}p/G  V^G  G°°(P/G). 

Here  /j([x])  =  P(x)  for  any  equivalence  class  [x]  in  P/G. 

From  the  definition  of  the  characteristic  distribution  on  a  Poisson  manifold,  it  follows 
that  the  hamiltonian  vector  field  Xh  leaves  invariant  the  symplectic  leaves  of  P/G.  Thus 
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any  Casimir  function  on  P/G  is  an  integral  of  motion  for  Xu-  The  integral  curves  of  A  u 
project  under  tv  to  integral  curves  of  X /, . 

The  process  sketched  out  so  far  is  the  essence  of  Poisson  reduction.  See  [9,  10]  for 
details  and  generalizations. 

For  our  purposes  the  main  examples  of  Poisson  manifolds  are  T*G  and  Q* .  The  Poisson 
structure  on  T*G  is  given  by  the  canconical  symplectic  structure  u>  on  T*G  as, 

{4>i  V’It’G  —  w  (X^,  X^)  (11) 

where  io  is  the  canonical  symplectic  form  on  T*G  (c.f.  Abraham-Maisclcn  p.  315). 

The  second  example  of  interest  to  us  is  the  space  Q*  the  linear  dual  of  a  Lie  algebra 
Q.  The  space  Q *  carries  the  Lie-Berezin-  Kirillov  -  Kostant  -  Souriau  Poisson  structure  (s), 
defined  by 


{<£,  i/’It  (h)  =  T  (h> 


8<j> 

8il> 


(12) 


where  (j>,  ij)  €  C°°(Q*)  and  /(  6  Q* .  Here  the  variational  derivative  is  defined  via 
the  Frechet  derivative: 


DM- 6  =  ,  64} 

for  6  <f>  €  Q* .  By  finite  dimensionality  ,&**  —  Q  and  €  Q. 

The  minus  (plus)  bracket  in  obtained  by  viewing  Q*  as  the  left  (right)  Poisson  reduction 
of  T*G  by  G,  as  explained  further  in  the  next  section. 

3.  Lie  -  Poisson  Reduction  and  the  Maximum  Principle 

Recall  that  the  tangent  and  cotangent  bundles  of  Lie  groups  are  trivial  as  vector 
bundles,  i.e. 


TG  ~  G  x  g  and  T*G  ~  G  x  Q*. 
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The  trivializations  are  made  explicit  by  tangents  of  left  and  right  translation  operators. 
Thus 


(id,  A)  :  T*G  ->  G  x  Q* 

ag  ^  (g,  TeL*ag ) 

(id,  p )  :  T*G  G  x  Q* 

ag  ^  (9,  TeR*g  ag ) 

are  the  corresponding  Poisson  diffeomorphisms.  In  the  case  of  G  =  50(3)  the  rota¬ 

tion  group,  these  diffeomorphisms  correspond  to  passage  to  body  and  space  coordinates 
respectively,  in  rigid  body  mechanics. 

Observe  that  the  canonical  symplectic  structure  on  T*G  is  G'-invariant  and  is  best  ex¬ 
pressed  in  one  of  the  trivializations  (A  &  M  p.  315).  In  each  case  the  quotient  T*G/G  ~  Q* . 

Proposition  (Lie-Poisson  Reduction): 

The  quotient  maps 


p  :  t*  G  —>  Q* 


are  Poisson  homomorphisms  if  for  A  we  take  the  minus  Lie-Poisson  structure  on  Q*  and  for 
p  we  take  the  plus  Lie  Poisson  structure  on  Q* . 

Proof.  It  is  a  calculation  using  the  expression  of  u  on  T*G  in  one  of  the  two  trivializations, 
and  using  the  definition  of  quotient  Poisson  structure  as  in  section  2.  See  for  instance  [9] 
[10]. 

Denote  the  space  Q*  with  the  minus  (plus)  Poisson  structure  (12)  as  Q*_  (Q\).  From 
the  discussion  in  section  2,  every  left  -  invariant  hamiltonian  on  T*G  gives  rise  to  a  reduced 
dynamics  on  Q*_.  (For  our  optimal  control  problem,  left  reduction  is  of  interest.) 

Suppose  h  :  Q*_  — >  1R  is  a  reduced  hamiltonian  function.  Let  {Xl,  •  •  • ,  A’„ }  be  a  basis 

for  the  Lie  algebra  Q.  Let  {Xj,  •  ■  • ,  }  be  the  dual  basis  for  Q* .  Thus  (A'?b,  Xj)  —  bl- 
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the  Kronecker  Symbol.  Any  p  £  Q*_  can  be  expressed  as  /<  =  pi  A”b.  The  differential 

(=1 

equations  satisfied  by  the  coordinate  functions  are 
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A«  —  {  i  i  h } 


*  =  1,2,-  •  •  ,n 


(13) 


where  the  Poisson  bracket  has  to  be  computed  from  (12).  Explicitly 


n  n  ci 

Ai  =  -  XIS  flk  Fb  2’  (14) 

j=i  fc=i  tij 

Here  rfj  are  the  structure  constants  of  Q  in  the  basis  {Ad,  •  •  • ,  Xn}. 

By  the  properties  of  the  Poisson  bracket  it  is  immediate  that  h  is  constant  along 
trajectories  of  (13).  Additionally,  the  Casimir  functions  on  Q*_  are  also  constant  along 
trajectories  of  (13).  For  Q  —  50(3)  the  Lie  algebra  of  skew  symmetric  matrices,  any 
Casimir  function  is  of  the  form  $  (//{  +  /A  +  fi\).  In  general  there  may  not  be  any 
nontrivial  (i.e.  non-constant)  Casimir  functions.  (The  Lie  algebra  of  the  affine  group  or 
ax  +  b  group  is  an  example.) 

Returning  to  the  optimal  control  problem  of  this  paper,  since  the  hamiltonian  H  in 
(9)  is  already  expressed  as  function  on  Q* ,  we  note  that  the  reduced  hamiltonian  is 

h  =  (f,  (o)  +  (  E  (15) 

l~  I 

We  have  in  effect  shown  the  following  reduction  of  the  maximum  principle. 

Theorem.  Consider  the  optimal  control  problem, 


subject  to 


Min 

u 


It  uj.dt 


g  —  TeLg  ■  I  £0  +  ^  ui  {t)£i 


i=i 


Sf(0)  =  go  and  g(T)  =  gx. 
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Then  every  regular  extremal  is  given  by 


Ui  =  (/x,  6)  /  C, 

where  fi  is  an  integral  curve  of  the  vector  field  Xh  on  Q*_  corresponding  to  the  hamiltonian 

h  =  fe.  M  + 

i=  1 

and  the  Poisson  bracket  {•,•}_  on  Qt_  is  given  by  (12),  In  coordinates  on  Q*_  the  integral 
curves  satisfy  the  ordinary  differential  equations  ( 13 )-( 14 ).  ■ 

4.  Example:  The  Unicycle. 

Consider  the  unicycle  equations 

X  —  COS  (</>)  U2 

y  =  sin  ((f))  u.2 
<j)  =  ui 

These  equations  are  equivalent  to  the  nonholonomic  constraint  x  sin  4>  —  ij  cos  0  =  0. 
Set, 
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COS  (f) 

sin  (f> 

0 


—sin  4> 


COS  4> 

0 


X 

y 

l 


clearly  g  G  SE( 2)  the  rigid  motion  group  of  the  plane.  Further  the  unicycle  equations 
take  the  form 


9  -  g  •  PC  u i  +  X 2  u2) 


(16) 


where, 
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0-10 
1  0  0 

0  0  0 


Xx 


x2 

clearly  this  is  a  left-invariant  system  in  SE(2).  Further,  since  the  Lie  algebra  of  SE(2)  =  $e(2) 
=  span  {Xi,  X2 ,  -X^3  } ,  where 

‘0  0  O' 

X3  =  001 

0  0  0 

and  [Xi,  X2]  =  X3,  (also  [Xi,  X3]  =  —  X2  and  [A'2,  A' 3]  =  0),  it  follows  that  we  have 
a  controllable  system  on  SE( 2).  (This  way  of  looking  at  the  unicycle  was  pointed  out  by 
Naomi  Leonard.) 

Next  consider  the  optimal  control  problem 

T 

min  J  (11%  +  u\)  dt 

0 

subject  to  end-point  conditions.  By  our  theorem,  regular  extremals  are  given  by  integral 
curves  of  a  reduced  hamiltonian  on  se(2)l.  The  hamiltonian  is  in  fact,  in  the  coordinates 
corresponding  to  the  dual  basis  {Aj,  Xij,  A3  } ,  given  by  h  —  (/<j  +  //J)/2.  The  Poisson 
bracket  of  two  functions  <f>  and  if)  is  given  by 

{</>,  V>)-  =  V^T  A  W  V  </, 

where, 

3 

[AM],-;  =  T  r«  w- 

fc=l 
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Explicitly, 


/  0  -n  3  H2  \ 

A  =  Ms  0  0  .  (18) 

\-/i2  o  o  J 

The  reduced  hamilton’s  equations  are  (from  (14)), 

Mi  =  -M2M3 
M2  =  MlM3 

M3  =  MiM2-  (19) 

The  Casimir  functions  are  of  the  form  =  $  (//J  +  // 3 ) ,  (equivalently  V$  is  in  the  kernal 
of  A).  The  level  sets  of  Casimir  functions  (i.e  symplectic  leaves  in  <3*_  )  are  right  circular- 
cylinders  {fi  :  P2  +  // 3  =  c}.  Integral  curves  of  (19)  are  intersections  of  level  sets  of  h, 
also  cylinders  { \i  :  =  2 h}:  with  the  symplectic  leaves.  Thus  the  phase  curves 

on  symplectic  leaves  are  as  in  Figure  1. 
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Equations  (19)  are  explicitly  solvable  by  elliptic  functions.  To  see  this,  note  that, 


fi2  +  (2 H  +  C)  /J2  -  2 /fil  =  0.  (20) 

This  is  the  equation  of  an  anharmonic  oscillator  with  quartic  potential  term.  The  general 
solution  to  (20)  is  given  by, 


fi 2  (t)  =  P  Sn  (A (t  -  t0),  k) 

where  Sn  (u,  k )  is  Jacobi’s  elliptic  sine  function,  A  s.t.  A2  <  (2 H  +  C)  <  2A2,  to  is 
arbitrary,  k 2  =  2H At  C  —  1,  ;32  =  2 H  +  C  —  A2.  Then  // 1  and  // 3  are  determined 
from  Hi  —  a/2 H  -  h\\  A4 3  =  \/C  -  //2 . 

The  optimal  controls  are  given  by 

«i  =  /<T,  U2  =  H2 ■ 

It  can  be  shown  that  in  this  example,  there  are  no  abnormal  extremals  since  the  strong 
bracket- generating  condition  holds  (c.f.  Strichartz  [12]) 

5.  Conclusion 

We  have  worked  out  explicitly  the  Poisson  reduction  of  certain  G'-invariant  optimal 
control  problems  on  Lie  groups.  The  approach  presented  here  yields  an  algorithm  for 
constructing  regular  extremals. 

6.  Notes 

This  paper  is  meant  to  be  used  as  a  set  of  notes  for  a.  mini-lecture  at  the  Workshop 
on  Mechanics,  Holonomy  and  Control,  IEEE  CDC  ,  Dec  14  1993.  The  main  ideas  of  this 
paper  are  already  to  be  found  in  the  papers  of  Faybusovich  referred  to  below.  He  used 
the  language  of  collective  hamiltonians  and  we  use  the  language  of  Poisson  reduction.  We 
are  also  a  bit  more  fussy  about  right  and  left  reduction  and  the  different  reduced  Poisson 

structures  one  associates  with  them.  In  the  papers  of  Jurdjevic,  the  reduction  is  carried  out 

in  various  examples  by  hand,  but  Jurdjevic  tends  to  work  with  the  symplectic  structures. 
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These  papers  of  Jurdjevic  and  the  references  therein  contain  many  additional  examples  of 
explicitly  solvable  optimal  control  problems.  Our  own  enthusiasm  for  Lie-Poisson  structures 
dates  back  at  least  to  our  work  [7]  on  the  dual  spin  problem.  This  and  later  collaborations 
have  influenced  our  thinking.  A  special  thanks  is  extended  to  Leonid  Faybusovich  for  useful 
and  enthusiastic  discussions. 
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